=3 

1—5 



INITIAL TIME SINGULARITIES IN NON-EQUILIBRIUM EVOLUTION OF 
CONDENSATES AND THEIR RESOLUTION IN THE LINEARIZED 

APPROXIMATION 

J. Baacke('^ '=), D. Boyanovsky(^''=) , H. J. de Vega(= ^) 

(a) Institut fiir Physik, Universitdt Dortmund, D-44221 Dortmund, GERMANY 
(b) Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260 USA 
(c) LPTHE, Universite Pierre et Marie Curie (Paris VI) et Denis Diderot (Paris VII), Tour 16, ler. etage, 4, Place Jussieu 

75252 Paris, Cedex 05, FRANCE* 
(February 1, 2008) 

The real time non-equilibrium evolution of condensates in field theory requires an initial value 
problem specifying an initial quantum state or density matrix. Arbitrary specifications of the initial 
quantum state (pure or mixed) results in initial time singularities which are not removed by the usual 
' renormalization counterterms. We study the initial time singularities in the linearized equation of 

On . motion for the scalar condensate in a renormalizable Yukawa theory in 3 + 1 dimensions. In this 

renormalizable theory the initial time singularities are enhanced. We present a consistent method 
for removing these initial time singularities by specifying initial states where the distribution of 
high energy quanta is determined by the initial conditions and the interaction effects. This is done 
through a Bogoliubov transformation which is consistently obtained in a perturbative expansion. 
The usual renormalization counterterms and the proper choice of the Bogoliubov coefficients lead to 
a singularity free evolution equation. We establish the relationship between the evolution equations 
in the linearized approximation and linear response theory. It is found that only a very specific 
form of the external source for linear response leads to a real time evolution equation which is 
^ ' singularity free. We focus on the evolution of spatially inhomogeneous scalar condensates by 

I implementing the initial state preparation via a Bogoliubov transformation up to one-loop. As a 

. concrete application, the evolution equation for an inhomogenous condensate is solved analytically 

fT^ ' and the results are carefully analyzed. Symmetry breaking by initial quantum states is discussed. 

o 

ON 

2 '' The study of the real time dynamics and the evolution of non-equilibrium quantum states has now become ubiquitous 
(~| in cosmology and high/intermediate energy physics. In cosmology the real time evolution of expectation values of 
IL" quantum fields is a necessary component of a microscopic description of the inflationary dynamics and the subsequent 
hot FRW stage (big bang) seeking to give a realistic description of the early universe and the physical processes 
originated there. In the physics of heavy ion collisions a very active program seeks to establish potential experimental 
5-H ' signatures from possible non-equilibrium stages of the evolution of the quark-gluon and chiral phase transitions 
iQ. In cosmology a program that incorporates consistently the non-equilibrium evolution of initial quantum states 
or density matrices of thermal or non-thermal origin including renormalization and backreaction effects had been 
pursued vigorously during the last few years In high/intermediate energy the possibility of studying the quark- 

gluon plasma and chiral phase transition at the forthcoming ultrarelativistic heavy ion colliders (RHIC and LHC) 
has motivated a substantial effort to study out of equilibrium dynamics during phase transitions. In particular the 
formation of coherent pion domains H, the evolution of non-equilibrium initial density matrices and states of high 
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I. INTRODUCTION 



energy density |9(] , and isospin condensates [10| . Non-perturbative techniques had been developed to study consistently 
non-equilibrium dynamics of quantum field theories and current computational facilities allow the possibility of 
studying the non-equilibrium dynamics of non-linear, inhomogeneous configurations in quantum field theories |Tl| ] 
including gauge theories p^JT^ , for which recent lattice simulations of non-equilibrium gauge field theories with 
topological excitations had recently been reported . 

The real time evolution of either density matrices or pure states, or alternatively of matrix elements must be set 
up as an initial value problem, either by specifying the initial state or by providing the Cauchy data (expectation 
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values and their derivatives) typically on space-like hypersurfaces. Once this initial value problem has been set 
up at some initial time, the real-time evolution of the expectation values or other matrix elements can be studied 
either analytically in the case of small amplitudes [ p"5|Jl^ or numerically in the case of large amplitude configurations 
|p|,P| jTl|jr^ , [l^ (Although analytic expressions are available in some extent). 

An important but largely unnoticed subtlety arises in these situations in that besides the usual ultraviolet diver- 
gences associated with masses, couplings and wave-function renormalizations there appear initial time singularities 
|^,|8|. The physical reason for these initial time singularities can be understood as follows: the initial state (either 
pure or mixed) is typically chosen to reflect some physical description but generally is either some initial pure excited 
state with free field quanta or a thermal density matrix for free field theory. The choice of the initial state (including 
the field expectation value and its time derivative) has been essentially arbitrary and in particular independent of 
the field hamiltonian. The time evolution with the interacting Hamiltonian suddenly couples at the initial time the 
infinite number of degrees of freedom of the theory, redistributing the spectral densities. In the case in which the 
underlying theory is renormalizable this redistribution of the spectral densities results in a divergent response. Such 
effect is also present on systems with a finite number of degrees of freedom but it is then finite. 

The consideration of singularities associated with setting up initial conditions in a quantum field theory has been 
addressed originally by Stueckelberg |^ , the similarity with sharp boundary conditions in a Euclidean formulation has 
been studied by Symanzik [|9| and has since found different possible solutions |p|,^| jl7|Jl^ , |20|] . A very appealing method 
to prepare initial states that lead to evolution equations without initial time singularities has been recently advocated 



8[ for self-consistent real-time evolution. This method consists in defining an initial state as a Bogoliubov 



transformation of the initial states in a free field theory. The Bogoliubov transformation is chosen to cancel the initial 
time singularities. The advantage of this method is that it is physically transparent and can be implemented both for 
small amplitude, i.e. the linearized problem, as well as for the large amplitude case which must be necessarily studied 
numerically. 

In this article we focus on studying these initial time singularities and their resolution via the method of a Bogoliubov 
transformed initial state in a renormalizable theory in the case of small amplitudes of the scalar condensate. This case 
allows to obtain the evolution equations in a linearized approximation with an analytical solution to the evolution. 
Furthermore, we establish the correspondence between this method and linear response theory for the case of linearized 
equations of motion of condensates. An important corollary of this correspondence is that only very specific choices 
of the external source in the linear response approach lead to a singularity free initial value problem. We choose to 
study the initial value problem for the evolution of a scalar field condensate in a Yukawa theory in 3 -I- 1 dimensions 
both for homogeneous as well as for inhomogeneous condensates. Whereas in fl^Jl^ the homogeneous case has been 
studied in a self-consistent manner and the linearized approximation has been extracted from it, the inhomogeneous 
case has not been studied, and hence we devote our attention mainly to this important case. In a renormalizable 
theory the initial time singularities are enhanced and new infinities associated with the preparation of the initial state 
emerge, this situation is highlighted in the renormalizable Yukawa theory which is the focus of our study. 

Main Results: i) The main results of our study can be summarized as follows: a definite proposal to set up the 
initial value problem in renormalizable quantum field theories based on an initially Bogoliubov transformed state. In 
order to eliminate the initial time singularities, the Bogoliubov coefficients are constrained to vanish in a particular 
manner for high momentum modes. In particular the Bogoliubov coefficients for a mode of momentum p must be 
chosen such that the 1/^,1/^"^ contributions to these coefficients are uniquely fixed by the initial data, the coupling 
and the mass [sec. IV]. Choices of Bogoliubov coefficients that differ by contributions of higher order in 1/p define 
different initial states, all of them free of initial time singularities. Thus, the time evolution of an initial state in 
quantum field theory is free of initial time singularities provided that high energy distribution of quanta of the initial 
state is specified in a very precise manner. 

This method is implemented consistently in the perturbative expansion and in combination with the usual renor- 
malization of mass, wave-function and coupling leads to a real-time evolution free of ultraviolet and initial time 
singularities, ii) As an example we study the real-time evolution of an inhomogeneous scalar condensate in the 
Yukawa theory, both in the case in which the scalar is heavy and can decay into fermion-antifermion pairs, and in the 
case in which the scalar is light and cannot decay into fermion pairs. Here we provide a detailed analysis of the real 
time evolution of an inhomogeneous scalar condensate corresponding to a spherical wave. 

The article is organized as follows; in section II we obtain the equations of motion for a scalar condensate. In section 
III we analyze the ultraviolet and initial time singularities. For simplicity we present first the case of a homogeneous 
scalar condensate. Section IV introduces the Bogoliubov transformed initial state in the case of a homogeneous 
condensate, discusses in detail the choice of the Bogoliubov coefficients that lead to an evolution free of initial time 
singularities and presents the singularity free real-time equations of motion for the homogeneous case. In section V 
we establish a relation between the initial value problem in the linearized approximation and linear response and 
discuss the constraints on the external sources that lead to a well defined initial value problem free of singularities. 
In section VI we extend the treatment to the case of inhomogeneous scalar condensates, obtain the corresponding 
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inhomogeneous Bogoliubov transformation consistently in perturbation theory and the equations of motion free of 
singularities to one-loop order. In section VII we obtain an analytic solution of the real-time equations of evolution 
for an inhomogeneous scalar condensate. We also provide a numerical analysis of the solution and discuss its main 
features. The conclusions summarize our work and discusses the potential applications of the methods presented. 

The Bogoliubov transformation of the tadpole and self-energy diagrams is presented in the Appendices A-C. Fur- 
thermore, Appendix D establishes several sum rules on the spectral densities. 



II. LINEARIZED EQUATIONS OF MOTION FOR CONDENSATES: 

Although the initial time singularities that will be discussed in this article are generic features of initial value prob- 
lems in field theory, they are highlighted in renormalizable theories. Therefore we choose to discuss these singularities 
and their resolution in a Yukawa theory in 3 -I- 1 dimensions. The focus of this article is to understand the problems 
of setting up an initial value problem to describe the non-equilibrium evolution of condensates or field expectation 
values in the linearized (small amplitude) approximation. Furthermore, we compare with an alternative formulation 
based on linear response. 

We consider a massive scalar field ^(x) coupled to a massive Dirac field ipix) in a Yukawa model specified by the 
Lagrangian density 

£($, ip, 4>) = ]^d^^{x) d^'^ix) + ^Af 2 -t- ij{x) [it^ + m + g <^{x)] ip{x) . (Ill) 

We study the time evolution of the expectation value of the scalar field via the real time generating functional in 
terms of a path integral defined on a contour in complex time (CTP) ||2l|,^. The effective Lagrangian that enters in 
the contour path integral is 

£e// , 

where the ± labels on the fields refer to the forward (+) and backward (— ) branches corresponding to the forward and 



backward time evolution of the initially prepared density matrix. We now follow the procedure of references [J_5 16 
and use the tadpole method to obtain the equation of motion for the expectation value of the scalar field 

(j){x) = < ^{x) > 

and write 

^^ix)^X^ix) + (t>ix) ; {x^ix))=0. 
We specify the initial data at the time io = by giving the initial condition, 

0(x, 0) = 0(x) and 0(x, 0) = ^{x) . 
Let us consider that the initial density matrix at time = is given by 

P(0)-|0>(0| 

with |0 > the free field Fock vacuum for the scalar and fermion fields. For t > 0, p{t) — e^*^* p{0) e*^* where H is 
the full Hamiltonian. This case is tantamount to considering an initial free field vacuum state and switching-on the 
interaction suddenly at t = 0. 

The equation of motion for 0(x, t) is obtained in a systematic perturbative expansion by imposing that (x, t)) — 
to all orders in perturbation theory. We will restrict our study to the case of small amplitudes of the condensate and 
will obtain the equations of motion linearized in 0(x, t) . In this linear approximation the self-energy kernel is obtained 
to any desired order in a perturbative expansion in the Yukawa coupling but in the state with vanishing condensate. 
Thus assuming that the state with vanishing condensate is spatially translational invariant it is convenient to perform 
a spatial Fourier transform for the condensate and the self-energy kernel |jl^,0 . Anticipating renormalization effects 
we introduce the renormalized field and mass in the Lagrangian before shifting the field by the condensate 

Since the fermionic fields will be integrated out to obtain the equation of motion for the expectation value of the 
scalar field, we do not introduce the renormalizations associated with the fermionic fields. We now drop the subscript 
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R from the renormalized quantities to avoid cluttering of notation, with the understanding that the scalar field and 
its mass are the renormalized ones. 

The equation of motion for the spatial Fourier transform of the condensate 



to one loop order is given by (see ref. [Tm for details) 



(1 + SZ) 4>q{t) + (t>q{t) + {M^ + SM^) <j)q{t) + / dt' Sq(t - t') (j)q{t') - J = 



(II.2) 



with SZ = Zj, — 1 and 



/(Pk 
trS>{t, t) 
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(2^ 



tr [S>{t-t')S<_^{t' ^t)-S<{t~t')S>_^{t' -t)] . 



(11.3) 
(11.4) 



The contribution J is due to the fermion tadpole, it is usually absorbed in a constant shift of the field (t){x). The 
fermionic Green's functions in the vacuum state are given by |16| 



S>{t,t') = -i / d3a;e-'P-^ (V'(x,i)V;(0,t')) 



2Er, 



S<{t,t') d^x e-'P-'' (V;(0,i')V'(x,t)) 



2E„ 



-^EAt^t') + ^) + e-*^''(*-*')7o(|( - m)7o 



Ep = y/p^ + vn? 



III. ULTRAVIOLET RENORMALIZATION AND INITIAL TIME SINGULARITIES 



The evolution equation of the type (II.2) contains two types of divergences: i) ultraviolet divergences which are 
removed by the mass and wave function renormalizations, ii) initial time singularities. 

To illustrate these singularities in a more clear manner we now focus on the case of homogeneous condensate, i.e. 
q = 0. We find 

I^o(t-t') ^-2,^/ sin[2i.,(.-0] , (III.1) 



J = -4 m 



{2TiY2Ep 



Whereas J acts as a constant source term and can be absorbed in a shift of the expectation value 4>{x \ th e self-energy 
kernel leads to ultraviolet divergences as can be seen upon integrating by parts the non-local term in (tL3) three times 



* 2 f d^p rt 







dt' S]o(i - t') Mt') = -g' J ^2n/2Ep Jo 'i^[2i?p(i - t')] Mt') = 

r (fir) f 1 1 1 • 

I J2^i; (2^ ^^^'^ ^ 2^ 0o(O)cos(2i.,t) - ^ ^o(O) sm2i.,. 

1 .. 1 .. 1 /■* 

Mi) + TTTTTTT '/'o(O) cos2Ept + / dt' cos[2Ep{t - t!)\ </>o {t') 



{2EpY ' {2Ep)^ ^ " ' {2EpY Jo 
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Using dimensional regularization the coefficient of (f)o{t) becomes 



{2Tr)^2Ep Ej 47r2 



1 , 47ru2"l 
7 + - + In —9 



This agrees with the expression obtained by evaluation the corresponding Feynman graph in 4 — e dimensions. The 
renormalization is performed here at = 0. 

The coefficient proportional to <f)o {t) is the wave function renormalization which again in dimensional regularization 
is given by 



6Z = i^ 



{2TT f2Ep £4 8^2 



- -7- 
e 



2 , Att^i^ 
^ + In 



These ultraviolet renormalizations arc cancelled by the mass and wave function countcrterms but there still remain 
singularities arising from the terms that are evaluated at the initial time t = Q, these are given by 



dt' ^o{t-t') Mt') 



f 

Jo 

Wj^i; ^^^^^^^ 

+ (2^ "^"^^^ cos2i;pi 

Obviously this expression is singular as t ^ 0. Simple power counting shows that the coefficients of 00 (0), '/'o(O), 
and 00 (0) diverge as 1/t, and logt, respectively. At finite t they are finite due to the oscillatory behavior of the 
integrand. 

The physical reason for these singularities is the following: having prepared the initial state to be the free field 
Fock vacuum and switching-on the interaction suddenly at the initial time t — 0, the interaction redistributes the 
spectral density of fields. The scalar field states overlap with the fermionic continuum of states and the particles 
become dressed by the interaction. This dressing effect which is responsible for mass, wave function and coupling 
renormalizations occurs suddenly when the interaction is switched on and is refiected as an initial time singularity. 
Obviously these short time singularities will be present at finite temperature or density, and are a consequence of the 
fact that the the underlying field theory posses an infinite number of degrees of freedom. 

Our main point in this article is that these initial time singularities can be removed and an initial value problem 
can be defined consistently and free of singularities by considering an appropriately chosen initial state that is dressed 
by the interactions. We thus propose to initialize the real time evolution by providing an initial state that includes 
the dressing as a Bogoliubov transformation from the free field Fock states. Furthermore we will argue that this 
construction leads to a consistent initial value problem for real-time dynamics and can be implemented systematically 
order by order in perturbation theory. 

We now discuss in detail this procedure in the example under consideration to lowest order in the Yukawa coupling. 



IV. EQUATIONS OF MOTION WITH BOGOLIUBOV TRANSFORMED STATES: HOMOGENEOUS 

CASE 

In this section we introduce the Bogoliubov transformed states and show explicitly how the introduction of these 
dressed states provides a solution to the problem of initial time singularities. For the sake of clarity we study first the 
homogeneous case q = and postpone to a later section the generalization to inhomogeneous condensates. 

From the free field Fock vacuum state |0) a Bogoliubov transformed state is obtained after a unitary transformation 

m = e-Q\o). 

with Q a unitary operator. Bogoliubov transformed operators are defined via 

Ob = cxp(-Q) O cxp(Q) , 

therefore if the vacuum state |0 > is annihilated by the destruction operators, the Bogoliubov transformed annihilation 
operators annihilate the state |0|,). 
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To lowest order in the Yukawa coupling the Bogoliubov transformation that required to cancel the initial time 
singularities only involve fermionic fields. Only when scalar contributions in higher order corrections arise there will 
be a need to introduce the Bogoliubov transformation for scalar fields. 

Therefore we introduce the antihermitian operator Q that generates the unitary Bogoliubov transformation 



This generator of Bogoliubov transformations illuminates at once the nature of the Bogoliubov transformed initial 
states. Acting on the free field Fock vacuum state, the Bogoliubov transformation leads to a state that is a linear 
combination of particle-antiparticle pairs ot total zero momentum. The fact that the total momentum of these pairs 
is zero is of course a result of the fact that the scalar condensate is homogeneous. 

The commutators of Q with the fermionic creation and annihilation operators are given by 

[Q,6(p,.s)] =0ps C^'"" dH-p,s) 
[Q,d\-p,s)\ =-/?y,e-*^-6(p,s). 

which leads to the following relation between the transformed and the original operators 

6(p, s) = cos/3p,s 6|,(p, s) + sin/3p,s e"'"'-' rfj(-p, s) , 
dt(-p,s) = -sin/3p,s e"'*"'^ 6^ (p, s) + cos /3p,s dl{-p,s) . (IV.l) 

The operators 6(p, s), rf(p, s) as well as 61, (p> s), d\,{Tp, s) obey the usual canonical anticommutation relations. 
The initial density matrix is now given by 

Pb(0) = e-« p(0) e« = |0b >< Obi . (IV.2) 

Although we have focused on a pure (vacuum) state, obviously this can be easily generalized to thermal or non- 
thermal mixed states. 

In Appendix A we provide the details that lead to the following Green's functions in the Bogoliubov transformed 
states. 

With this restriction the transformed Green function becomes 

i5>(t,x;t',x') = (Ob|V(t,x)Vi(i',x')|Ot) = 



/ 



d^p 



and 



-sin/3p cosppe''^" ^p(j/ + m)j5jae-'^''^*+*'^ 
- sin/3p cos l3p e"'^" S p 7570 m) e*^''(*+*') 
+ sin2/3p757o(^+m)757o e*^''^*"*')] , 



i5<(i,x;t',x') = (Ob|Vi(t',x')V(i,x)|Ob) = (IV.3) 

d^p 



sin/3p cos/3p e'^^^ Sp^+ 771)7570 e-'*^''(*+*') 
+ sin/3p cos/3p e-'^^' Sp 7570 m) e'^^^*+*'^ 



-cos^ 



/3p 757ofe/+m)757oe*^''(*~*') 



Perhaps the most striking feature of these Green's functions is their lack of time translational invariance, the main 
reason is that the Bogoliubov transformed states are not eigenstates of the bare particle number. We note that in 
the terms with t + f' 'a translation of the time variables can be compensated by a change in the phase 5p, i.e. a gauge 
transformation of the fermionic fields. 
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It is this lack of time translational invariance that will allow to cancel the initial time singularities as shown explicitly 
below. 

Th e evolution equation is obtained in the same manner as in the previous section and is exactly of the same form as 
(II. 2) with q = as befits the homogeneous case, with the self-energy and tadpole kernels now given by the expressions 
(II.3)-(II.4) but in terms of the Bogoliubov transformed Green's functions 



and 



Ji,{t) = -ig tr S'^>(t,x;t,x) = -9 J (^2-Kf2E t"^"" 2/3p - {-Ap) sin2/3p co&{2Ept ~ 5p)] 

Eb,o(i, t') = ~2g^ j (2^)32i^^ ^ {« cos(2/3p) suv[2Ep{t - t')] 
-8 p msin(2/3p) [sin(2£'pt - 5p) - s\n{2Ept' - 5p)]} . 



(IV.4) 



(IV.5) 



The first, time independent contribution to J\){t) in cq.(IV.4) can be absorbed into a constant shift of the condensate 
much in the same way as in the un-transformed case. The second, time dependent term will be used to cancel the 
initial time singularities. 

The self-energy kernel ( IV.S ) can be written as So(i — t') + AE|, o(i,t') with So(t — t') given by ( 111.1 ) and 

d^p 1 



AI]b,o(t,i')--2.9 



{8p^ (cos2/3p - 1) sin[2£'p(t - t') 



{2'KY2Ep 2Ep 
—8pm sin2/?p [sin(2i?pt — Sp) — sui{2Ept' — Sp)]} . 

We will assume in the following that (3p decreases sufficiently fast with p so that terms proportional to sin2/3p and 
cos(2/?p) — 1 lead to convergent integrals. This in fact will b e checked a posteriori when we find the required expression 
for Pp below. Then the ultraviolet divergences in eq.([V.5) are the same as in the perturbative vacuum. 

Integrating by parts now three times in t' in order to single out the ultraviolet and equal-time singularities from 
'^\>,o{tjt') in the equation of motion, we find that just as in the un-transformed case the ultraviolet divergences 
proportional to 0o(i) and (jjoit) are cancelled by the mass and wave function renormalization counterterms (up to 
finite parts depending on the renormalization prescription). 

The terms that result in initial-time singularities arise from 



JQ 

d^p 8p2 



(27r)32£;p Ep { 2Ei 



6o(0) cos2Ept 



{2Epf 



(/«o(0) sin2£'pt 



{2Ep) 



ho{Q) cos2Ept )■ cos2/3p 



We now require that these terms are cancelled by the time dependent terms of Jt,{t), eq.(IV.4). This requirement 
leads to the equation 



tan2/3p cos{2Ept — Sp) 



2p 

Ev 



2E„ 



ho{0) cos2Ept 



{2EpY 



/)o(0) sin2£'pt 



{2Epy 



^o(O) cos2Ept 



Comparing the terms proportional to the sine and cosine we find two equations that determine (5p , (3p. These 
equations can be solved perturbatively with 

Pp = bi,p 9 + b2,p + 0{g^) . 
To one loop order we only need to keep the linear term in g leading to 
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I3p cos Sp = g 



(ip sin 5p = g 



En 



E„ 



2Er, 



MO) 



{2Epf 



^0(0) 



{2Epy 



00 (0) 



(IV.6) 



At this stage we recognize that there is freedom in choosing the Bogohubov parameters to cancel the initial time 
singularities. The initial value problem will be free of singularities by choosing the coefficients Pp ; 5p so as to cancel 
the terms proportional to 1/p ,\/p^ and so that /3p vanishes as p — > cxd. 

Different choices of the Bogoliubov parameters that differ only in higher inverse powers of the momenta lead to 
different initial quantum states, but the initial value problem is free of initial singularities. This freedom is similar to 
choosing renormalization counterterms including finite parts, i.e. different renormalization prescriptions. 

Now the Bogoliubov correction and the mass and wave-function renormalization counter terms remove exactly 
all ultraviolet and initial time divergences from the equation of motion. To lowest order we can set /3p = in the 
self-energy in the equation of motion and having removed all ultraviolet and initial time singularities and absorbing 
the time independent contribution from the tadpole term in a constant shift of the condensate, we finally obtain the 
evolution equation in the case of the homogeneous condensate 



00 (t) + M' 00 (0 + / dt' ll,{t - t') 00 {t') = 
Jo 



where Yis{t — t') is the subtracted self energy kernel 



^s{t-t') = -g^ J 



d^p 



{2^Y2Ep E^ 



cos [2Ep{t-t')] 



(IV. 7) 



(IV.. 



The equation of motion above should be compared to Eqs (7.21) and (7.22) in [Q. The equation of motion obtained 
there differs from (IV. 7) by a finite renormalization that corresponds to a different renormalization scheme but 
otherwise the equations are the same up to one loop order. 

We note that even in the massless limit the tadpole J\, — (ipip) is nonvanishing, signaling chiral symmetry breaking, 
which in this case is solely a consequence of the initial conditions. This is an important result of our analysis, that 
non-equilibrium initial states can lead to symmetry breaking. 



V. LINEAR RESPONSE THEORY 



The initial value problem can be obtained by establishing direct contact with linear response theory as presented 
in ref. |2^ in the case of scalar condensates and in ref. |Q for fermion ic co herent states. This is achieved by coupling 
an external source term to the scalar field in the Lagrangian density (ll.l) 



V, ^] ^ C[<i>, V', ^] + Jextix) $(x) . 

The expectation value of the scalar field induced by this source term is given by 

($(x)) = I J dx' Jextix') [(<f+(x)$+(a;')) - {<f+{x)<f-{x'))] (V.l) 

where the superscripts ± refer to the forward and backward time branches in the real time generating functional. The 
bracket in (V.l) is the retarded commutator. As discussed in detail in ref. the inversion of (V.l) gives rise to the 
equation of motion for the scalar field with an inhomogeneity given by the external source term. Following the same 
steps detailed in the first section, we find the equation of motion for an homogenous condensate to be given by 



(1 + SZ) 00 (i) + (Af2 + SM^) Mt) + / dt' So(t - t') 0o(t') - J - Jext,o{t) 



which differs from (II. 2) in the lower limit in the non-local term with the self-energy. Assuming adiabatic s witch ing- 
on of the interaction from t = —oo we can now rewrite this equation of motion in a form that is closer to (11.2) by 
integrating by parts the non-local term. Definining, 

^"^^ " ^'^ = = :^^2,o(i - 1') = ^E^t - 1') 



dt'^' 



dt'^' 
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and integrating by parts we obtain 

(1 + SZ) Mt) + {M^ + 5M^) Mt) + Si.o(O) 0o(<) - S2,o(0) 0o(i) + ^3,0(0) Mt) 



I dt' E,,o(i - t') 4>o {t') - J = JextAt) : 

^ — 00 



with Ss.o(^ — t') given by (IV.S). Using the exphcit form of Xlo(t) given by ([II. 1) we find 



-6ivr 



S2,o(0) = 



= -5Z 



Therefore the specific choice of external current 



Jext{t) 



dt' ^s{t-t') '4>0 {t') 



(V.2) 



(V.3) 



(V.4) 



leads to the initial value problem described by eq.( |lV.7| ). 

This current depends on the past history of the condensate, and in general does not vanish for t > as it would 
be desirable from the point of view of linear response. In linear response the initial value problem is envisaged to be 
prepared by switching-on an external source and the interaction adiabatically from t = —oo. The external source acts 
as a Lagrange multiplier, slowly displacing the condensate to the value to be determined at i = 0, and switching-off 
the source suddenly at this time. This allows the condensate to be formed and dressed over a very long period of 
time. The dressed condensate is then released when the external current is switched-off. However, in an interacting 
renormalizable theory this instantaneous switching-off of the external current results in singularities. To set up a 
consistent, singularity free and renormalized initial value problem as is the goal of this article, the choice of the 



current (V^4) is the one that establishes contact with a linear response formulation, in this case the current depends 
on the past history of the condensate, which obviously need not be specified for an initial value problem. On the 
other hand, such a choice of current, depending on the past history is rather artificial from the linear response point 
of view. 

We note that the current can be made to vanish at all times with the particular choice 

Mi) = 00 (0) + M0)t + \ MO) 

for t < 0. Obviously this behaviour manifests the problem of the initial time singularities as singularities in the 
behavior of the field at a remote past. 

An alternative would be to assume that the external source and therefore the condensate is adiabatically switched- 
on with a damping factor e*^* for t < Q but this results in discontinuities in the first or second derivatives aX t — Q and 
that would produce additional contributions from the integration by parts from these discontinuities. 

There are two main conclusions of this discussion on the relationship with linear response: 

• We have established a direct relationship between the evolution equations in the linearized approximation and 
linear response theory. The initial value problem free of UV and initial time singularities is shown to be 
obtained in the context of linear response through a particular choice of the external current. Such external 
current depends on the past history of the condensate and only a very specific form for it leads to a singularity 
free initial value problem. From the perspective of an initial value problem in which Cauchy data are specified at 
some given initial time on a space-like hypersurface this is a consistent choice. However, from the point of view 
of linear response this choice is somewhat artificial. The resulting external current does not vanish for t > Q. 
If we instead require an instantaneous switching-off of the external current at i = 0, initial time singularities 
become unavoidable. 

• The method of preparing a dressed initial state via a Bogoliubov transformation leads to a satisfactory description 
of the initial value problem. The usual mass, wave function and coupling constant renormalization counterterms 
cancel the ultraviolet divergences, and the Bogoliubov coefficients are judiciously chosen to cancel the initial time 
divergences consistently in perturbati on th eory. To lowest order in the Yukawa coupling and for a homogenous 
condensate such a choice is given by ( IV. 6 ). 

Having studied in detail the simpler case of the homogeneous condensate, we now move on to our main point, the 
study of the evolution of inhomogeneous condensates. 
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VI. EQUATIONS OF MOTION FOR INHOMOGENOUS CONDENSATES 



The equation of motion for non-homogeneous condensates in the amphtude approximation and in terms of spatial 
Fourier transforms reads 



(1 + 5Z) [^iq(t) + q2 <^q(t)J + (M2 + SW) Cj)^{t) 

+ / dt' Sq(t, t') (j)q{t') + 4,qW = . 

From the discussions of the previous sections we have learned that the Bogoliubov coefficients that define the dressed 
states at the initial time can be found consistently in perturbation theory. Since the self-energy is already of second 
order in the Yukawa coupling we will not need to consider the Bogoliubov corrections to the self-energy, but only to 
the tadpole term J,.q(f). 

The one-loop self-energy is therefore the usual one and given by 



{2TTf2Ep2Ep_^ 
(-2«)sin[(Sp + £;p_q)(t-t')] 



m 



(VI.1) 



The derivation of the tadpole diagram that determines Jti,q(0 using the Bogoliubov-transformed Green functions, 
is given in Appendix B. We find 



[ d'p 
J (27r)32^p 



d^p 



{2Tr)^2Ep2Ep-^ 



(q,p) e-'^^^+'^^-'>^* + j{(i,p') e'(-^p+-Sp-o) 



The function 7(q, p) is a function related to the angles of the Bogoliubov transformation, that will be specified below 
such as to remove the initial time singularities. 

The analysis of the singular and divergent contributions in the equation of motion proceeds as in the homogenous 
case, performing three integrations by parts with respect to the time in the non-local term, we find 



/' 

Jo 



dt' sin[(£;p + Ep_^)it - t')] 4>^{t') = 



MO) 



Ep + Ep-q Ep + Ep-q 



cos[(£;p + Ep_q)t] 



.(0) 



+ 



{Ep + Ep.qf 
{Ep + Ep_qf 



sm[{Ep + Ep-q)t\ 



lit) 



{Ep + Ep.^y 



cos[(£;p + Ep_q)t] + J ^ / dt' cos[(£;p + Ep_q){t - t')] <^q (f) . 



The parts containing 0q(O) and its derivatives lead to initial time singularities in the equation of motion, these can 
be isolated by writing 



/' 

Jo 



dt'Eq(i - t')<^q{t') 



Sing 



d^p 

{2TTf2Ep2Ep^q 



[EpEp_q + p(p - q) 



m 



r(q,p) e-^f^p+^p-o)* + r*(q,p) e'(^p+^p-'')* 



with 



T(q,p) 



.^q(O) 



+ 1 



.^q(O) 



+ 



9iq(0) 



Ep + -Ep-q (Ep -\- Ep-q)^ [Ep + Ep-q) 



note that (^q(O) 



(0). With the choice 
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7(q, P) = [Ep + Ep_q + p(p - q) - m^] T*(q, p) 

the current J\,^q{t) exactly cancels the initial time singularities in the non-local term with the self-energy. As in the 
homogeneous case, the current J|,,q(t), which is a Fourier transform of {'ip{x)tlj{x')) is nonvanishing. Again, chiral 
symmetry is here broken by the initial conditions. 

The ultraviolet divergent contributions of the self energy to the equation of motion are given by 



dt' Eq(t-0 </.q(0 



(27r)32£p2£p_q 

Mt) kit) 



UV div 

[Ep + Ep_q + p(p - q) - m'] 



Ep + Ep_q {Ep + Ep_q)3 
= Si(q2),^q(t) + S3(q2)?iq(i). 
Here we define the UV divergent parts of the self energy kernel in dimensional regularization as 



Ei(q2) = -8 5' 



S3(q')=8 5' 



d^-^p Ep + Ep_q + p(p - q) - 



(27r) 



3-e 



AEpEp-ci{Ep + Ep-d) 



d^-'p Ep + Ep^q + p{p-q) 



(27r)3 



AEpEp^qiEp + Ep^qf 



These expressions have to be regularized to obtain the renormalized equation of motion. This is discussed in detail in 
Appendix C. The singular and ultraviolet divergent parts are cancelled by the appropriate choice of the mass and wave 
function renormalization and the Bogoliubov coefficient in the tadpole. The final form of the subtracted self-energy 
kernel is given by 



S,,q(i - t') 



d^p 



{2TT)^2Ep2Ep. 



Ep + £'p_q + p(p - q) - 
{Ep + Ep.qf 



cos[{Ep + Ep_q){t - t')] , 



It follows from rotation invariance that the kernel 5]s_q(t) only depends on q^. 
The equation of motion in momentum space becomes 

[1 + 5Z + S3(q2)],^q(t) + + SZ) + M2 + SM^ + Ei(q2 



^{t) 



+ / dt' S,,q(t ~ t') <^q (t') = . 

Jo 

We show in Appendix C that we can decompose Ei(q^) and T,s{q^) as 

i;i(q2) = -SM^ -q^dZ + ASi(q2) 
X;3(q') = -<5^ + AE3(q2) 

The divergent and finite parts are explicitly given in Appendix C. We finally obtain the renormalized equation of 
motion which is free from ultraviolet and initial time singularities 



[l + AE3(q2)].^q(f)+ q2 + M2 + ASi(q2) cj>q{t) , 



+ 



f dt' 
Jo 



S.,q(i - t') 0q it') = , 



where 



E.,q(t-i') = -8fl' 



d^p 



EpEp-q + p(p - q) - 



(27r)34£;p£;p_q 



{Ep 



i?p-q)3 



COS [{Ep + Ep,)T] . 



(VI.2) 



(VI.3) 
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VII. SOLUTION OF THE EQUATION OF MOTION: NUMERICAL ANALYSIS 



We have derived in the previous section the renormaUzed equation of motion. It can be solved in a standard way 
via Laplace transform. We introduce 



V'q(s) 



dt e-'* (/.q(t) 



for the condensate, and 

for the self energy^. We find 



d^p 



d^p EpEp^d + P(P - q) ~ 

{2n)HEpEp^^ (Ep + Ep^^)^ 

EpEp^q + p(p - q) - m2 



cos [(£;p + Ep^q)T] 



{2Tr)HEpEp^^ {Ep + Ep^^f[{Ep + Ep^^f + s2] 
The Laplace transformed renormalized equation of motion becomes 



(VILl) 



[s" [l + AE3(q')+a,(s^q')] + q^ + + AEi (q^)} 



0q(O) + S(/.q(O) l + AS3(q') + ^s(s2,q2) 



(^s{s ,q j , 



so that 



0q(O) + 


50q(O) 




l + AE3(q2)+a 


s(s^q)' 


+ 0(O,q) a,(s2,q2)/s 


S2 


"l + AE3(q2)+a,(s2,q2y 


+ qM 


-M2 + ASi(q2) 



^q(s) 



The solution 0q(i) then is obtained by the inverse transformation 



zoo+c 



— 'tOC+C 



ds 



This solution is discussed in detail in Appendix D. The integral above is along the Bromwich contour with c a positive 
real constant to the right of all the singularities of the Laplace transform. The result can be written as (see Appendix 
D for details). 



0+ 



duj [cos(wt) 0q(O) uj Im Fi (— + io, q^) 



+ sin(wt) 0q(O) ImFi(-tj2 + io, q^) 
cosM ImF2(-c^2^io,q2) 



UJ 



where 



Fi(s2,q2) 
i^2(s^q2) 





1 - 


t-AE3(q')- 




q2) 


S2 


l + A^;3(q2)^ 


- CTs(s2,q2) 

|7s(s^q 


+ q2 4 


-M2 + AEi(s2,q2) 


S2 


l + A^;3(q2)^ 


- CTs(s2,q2) 


+ q2 4 


-M2 + AEi(q2) ' 



^The extra factor s is introduced so as to make as a function of . 
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0q(O) is not an independent initial value, it is determined by setting i = in the equation of motion (VI. 2) with the 
result 



<?iq(0) = -M^) 



q^ + M^ + AS]i(q^) 
1 + Al]3(q)2 



It is then convenient to define a kernel which combines Fi and F2 with prefactors uniqueley determined by 0q(O). 
With the definition 



1 q2 + AP + ASi 1 



q^ + + ASi 



the solution can be written as 

2 r 

<^q(0 = ~ / (?!)q(0) cos(wt) ImF3(cj — io, q^) + 0q(O) sin(ci;t) ImFi(— + io, q^ 

The function _F3 which is, up to prefactors, the Laplace transform of the solution, exhibits a pole at 

LOR ^ UJq + 5UJ ; UJq^ \/ A'P + q^ , 

with Slli = 0{g'^). If M > 2m the scalar field can decay into a fermion-antifermion pair, the pole actually describes a 
resonance. In perturbation theory the width of this resonance is perturbatively small and near the resonance we can 
approximate the function F3 by a Breit-Wigner resonance 



1 



Zk 



The resonance position is determined by 

ti;R = Re 

the residue is given by 



ujI + AEi 



1 + AS3 + aa(-w2 +io^q2) ' 



and the width by 



(l + AI]3)t^q d 



Fr = Im 

2a;R 



Rc 



< + ASi 



l + AI]3 + (5-,(-tj^ + io,q^ 



< + ASi 



1 + AE3 + a^(-w2_^io,q2) 



Numerical Analysis: 



We are now in conditions to study the evolution of an initial scalar condensate numerically by performing the 
inverse Laplace and Fourier transforms since all the quantities are given by the subtracted one-loop self-energy. 

We consider two separate cases: M > 2m in which case the scalar can decay into fermion-antifermion pairs, and 
M < 2m in which case the scalar particle is stable. In both cases we studied the evolution for an initial spherical 
wave of gaussian profile 

<^(0,x) = iVoexp(-xV2i?2) ; 0(O,x) = O ; with j <j)[id,^) = 1 (VII.2) 

Since this gaussian wave-packet has zero center of mass momentum, the peak of the wave packet will not displace 
under time evolution, but the wave packet will disperse and spread out in space-time. 

M > 2m: We have chosen M = 3m but there is a rather smooth variation of the wave function renormalization 
constant, position and width of the resonance for reasonable values of the ratio 2 < M/m < 10. As a first step we 
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calculated the value of the wave function renormahzation numericahy and find that Z^. differs from unity by less than 
3% and that the ratio of the width to the position of the pole F/wr « 0.02 for g = 1 or smaller and 2 < M/m < 10. 
Therefore if the coupling is of this order or smaller, the approximation of the spectral density (discontinuity across the 
cut) by the imaginary part of the Breit-Wigncr form given above is excellent. The agreement obtained from the full 
numerical evolution and that obtained from the approximate Breit-Wigner form is excellent. The evolution in this 
case is depicted in Fig. 1 that displays the profile of the condensate as a function of t and \x\ for A/ = 3; to = 1; g = 1 
. We see that the propagation is inside the light cone and damped in time. The peak at the origin is the center of 
mass of the wave packet, it oscillates in time with a frequency « and decays on a time scale F~^. The evolution 
at very early times is smooth. 

M < 2m: In this case the scalar field is stable, the spectral density features poles at w = iwR, below the fermion- 
antifermion cut. We find again that Zp « 0.97, implying, that the contribution from the fermion-antifermion contin- 
uum is negligible except for small t. Fig. 2 displays the time evolution of the gaussian wave packet in this case. The 
wave packet spreads in space-time, the evolution is always below the light cone as clearly illustrated in the figure. 
The amplitude of the wave packet decreases as a result of spreading. Its spatial integral, which is equal to (f>o{t), 
asymptotically oscillates with the pole frequency and amplitude Zp. 



VIII. SCALAR THEORIES 



The connection between the preparation of the initial state via Bogoliubov transformations and the formulation in 
terms of linear response allows to generalize the study presented above to scalar theories. In particular let us consider 
the case of a scalar self-coupled A^"* model using the linear response analysis. Focusing on the evolution equation for 
a homogeneous condensate to order we find 

(1 + 6Z) Mt) + (M' + SM^) Mt) + Si,o(0) Mt) + S3,o(0) Mt) 

+ [ dt' S,,o(t - t') 00 {t') = Je,t,o(0 . (Vlll.l) 



The order A tadpole has been absorbed into SM'^ as a mass renormahzation and th e self e nerg y to order A^ is given by 
the sunset diagram. The choice of counterterms 6Z, SM"^ are as given in equations ( V.2)-( V.3) and the external source 



introduced for linear response is specifically chosen as in eq.(V.4). In this manner we obtain the initial value problem 
free of ultraviolet and initial time singularities. An alternative interpretation of the external source associated with 
the linear response formulation can be provided by introducing a Hartree factorization in the Lagrangian with a term 
of the form 4 A $ {^^)- This term acts now as an explicit source in the linearized equation of motion which is chosen 
so as to lead to a singularity free initial value problem, thus requiring a non-trivial Bogoliubov vacuum. The linear 
response analysis leads very simply to a well defined initial value problem for a proper choice of the external source. 
The equivalence between the preparation of the state via a Bogoliubov transformation of the free field Fock vacuum 
(or density matrix) and the initial value problem obtained from linear response for the proper choice of external source 
allows now to generalize the results obtained above for the linearized approximation. 



IX. CONCLUSIONS 



The non-equilibrium evolution of condensates in real time requires to provide Cauchy data at some initial time, 
hence an initial value problem which requires the specification of an initially prepared quantum state or density 
matrix. An initial pure or mixed state of free field Fock quanta leads to initial time singularities. These have a simple 
interpretation: the evolution of expectation values or matrix elements in the interacting theory implies that the 
interaction is switched-on suddenly. The interaction rearranges the spectral densities of the fields and the response to 
the sudden switching-on of the interaction results in initial time singularities which are enhanced in a renormalizable 
theory. For ystems with a finite number of degrees of freedom such effects are also present but no singularities arise. 

In summary, the time evolution of arbitrary quantum states or density matrices in (interacting) field theory leads 
to short time divergences. Only for appropriately prepared pure or mixed initial states, as those considered in this 
paper, the time evolution is well defined. By appropriately prepared we mean states where the filling for high energy 
quanta follows a precise law determined by the initial data, couplings and masses. 

We have chosen to study these initial time singularities and provide a consistent resolution in a Yukawa theory in 
3-f 1 dimensions, this theory being renormalizable allows to identify all of the divergences and singularities: ultraviolet 
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divergences associated with mass, coupling and wave-function renormalizations and initial time singularities that 
cannot be cancelled by the usual counterterms. 

After recognizing the initial time singularities and their physical significance in the case of homogeneous conden- 
sates, we have proposed a rather simple approach to provide a singularity free initial value problem. We introduced 
Bogoliubov transformed initial states that incorporate the effects of dressing of states by the interaction. The Bogoli- 
ubov coefficients can be obtained in a systematic series expansion in the Yukawa coupling and we have obtained them 
to one-loop order in this theory. The usual renormalization counterterms cancel the ultraviolet divergences associated 
with mass, coupling and wave-function, and the Bogoliubov coefficients are chosen consistently to cancel the initial 
time singularities. That is, their high energy behaviour is fixed according to the initial data [sec. IV]. 

We have established contact with linear response theory by obtaining the evolution equations for the scalar con- 
densate and the initial value problem in the linearized approximation as the linear response to an external source 
coupled to the scalar condensate. This equivalent formulation clarifies at once the relationship between the linearized 
approximation for the evolution equations of the condensate and linear response. The corresponding initial value 
problem, i.e. providing Cauchy data for the field and its first derivative on a spatial hypersurface requires that the 
external source that couples to the scalar field does not vanish after the initial time. A very specific source term that 
depends on a given past history of the condensate furnished a singularity free initial value problem. 

After presenting the method in the simpler homogeneous case and establishing the relationship to linear response 
theory we focused on the important case of inhomogeneous condensates. Following on the steps for the homoge- 
neous case we have constructed the proper Bogoliubov transformation to lowest order in the Yukawa coupling and 
shown explicitly how a judicious choice of the Bogoliubov coefficients in combination with the usual renormalization 
counterterms leads to an initial value problem free of ultraviolet and initial time singularities. As an example of this 
consistent procedure we have provided a numerical study of the space-time evolution of an inhomogeneous scalar 
condensate both in the case in which the scalar can decay into fermion-antifermion pairs and in the case in which the 
scalar is light and stable. 

A noteworthy result of our study is the breakdown of (discrete) chiral symmetry via the initial conditions in the case 
of massless fermions. The initial state with a non-zero value of the condensate of its time derivatives (first or second) 
break the discrete chiral symmetry, we emphasize that this breakdown is not a consequence of an explicitly symmetry 
breaking term but of the initial quantum state with a non-equilibrium condensate. This is clearly a manifestation of 
non-equilibrium effects. 

We have also generalized the results of the fermionic case to scalar field theories by exploiting the relation to linear 
response, thus providing a generalized and consistent manner of describing non-equilibrium evolution of condensates 
in terms of an initial value problem free of ultraviolet divergences and initial time singularities. 

Applications: We foresee several applications of these methods: i) we can now study consistently the evolution 
of inhomogeneous pion condensates after a chiral phase transition by setting up a physically reasonable initial value 
problem that incorporates the important features of the transition in the fully interacting initial state (or density 
matrix). This approach is complementary to that advocated in ref. ii) In cosmology we can now study the non- 
equilibrium dynamics of inhomogeneous configurations by providing the initial field profile and the first derivative on 
a space-like hypersurface and following the space-time evolution of this configuration. In particular a very relevant 
setting for cosmology is that of supersymmetric theories during for example the stages of rolling of the scalar field 
component. Treating the dynamics as an initial value problem, the initial conditions on the scalar field, displaced 
from the equilibrium position breaks supersymmetry. This breakdown of supersymmetry is not explicit at the level of 
the Lagrangian, but by the quantum state. Our formulation allows us to follow the dynamics consistently and study 
the consequences of this supersymmetry breaking. Work on these issues is in progress. 

The next step in our program is to extend the results obtained in this article, valid in the linearized approximation, 
to a full non-linear inhomogeneous problem. We expect to report on progress on these and other issues in the near 
future. 
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APPENDIX A: FERMIONIC GREEN'S FUNCTIONS IN THE BOGOLIUBOV STATE FOR THE 

HOMOGENEOUS CASE 



The Green function ^-^(t, x; i', x') in the Bogoliubov-transformed state jO^) is defined via 



i5>(t,x;i',x') = (Ob|V(t,x)V;(t',x')|Ob) 
and the (free) field operators are given, in terms of the creation and annihilation operators, by 



(A.l) 



{2'KY2Ep, 



6t(p',s') [7(p',s') e- 



-iE.t' 



di-p',s') V{-p',s') e^^^'*' 



(A.2) 



We use the n orma liz ation fo r the spino rs U jp, s)U{p, s) — 2m so that W{p,s)U{p,s) = 1. 

Using eqs.( [V.l ), ( 1V.2| ), (A.l) and ( [A.2| ) we obtain the following expression for the transformed Green function 
5>(i,x;i',x'), 

iS>{t,^;t',^') = (Ot|V(i,x)V;(t',x')|Ob) = 



E 



{2ny2Ep 



„ip(x-x') 



'cos\(3ps)U{p,s)U{p,s) e-^^-(*-*') . 



-sin/3p, cos/3p, e'*- U{p, s)V{-p, s) e-'^''(*+*') 
-sin/3p, cos/3p, e-^*-U(-p,s)i7(p,s) e*^''(*+*') 
+ sin^ F(-p, s) V{-p, s) e'^-^*-*') . 

As long as we consider homogeneous condensates, the angles Pps and 6ps can be chosen to depend only on the 
modulus IpI, and on the helicities. The weight of the two possible helicities is still arbitrary and we consider these 
angles to be functions of the helicity matrix 



Sp = 



crp 
crp 



where p = p/|p|. We have, e.g., 

/(Sp) E U{p, s)V{-p, s) = E t^(P' «)^(-P' «)/(^P) = E fi^)U{P, s)V{-p, s) . 

s s s 

Sp commutes with the all other matrices that are relevant to this discussion, 70,75 and P7 and can therefore be 
treated as a c-number. 

Specifying the vector p'' to be on shell, p^^ = (Ep, p), and using 70770 — —7 it is straightforward to find 

J2uip,s)Uip,s) =]J+m, 

S 

E^(-P;'5)t^(-p, s) ==7o(x/-™)7o = 757o(x/+ »77.)757o • 



the mixed product can be found by resorting to the representation 



U{p,s) = 



^Ep+7 



V{-p,s) 



10(1/ - '7i)7o / 



Ep + m \ X 
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Where Xs is an eigenspinor of crp with eigenvalue s from which we find 



X! ^(~P' s)U{i), s) = 7570 m) 

S 

U (p, s)V{-p, s) = m)757o 



We note that the phase 6ps appears in a combination such that a shift in this phase can be compensated by a shift 



in the origin of time, i.e. a time translation t t + to Sps 



2Epto. Since the square of the helicity matrix 



is the identity and the only odd function of Pps multiplies the mixed terms, we found that the simplest Bogoliubov 
transformation that is required to cancel the initial time singularities is such that the phase Sps is independent of s 
and that f3ps = Sp/3p so that 

cos f3ps = cos Pp 
sin /3ps = Spsin/3p . 

Such a choice has proven to be appropriate for removing the initial singularity for the full one-loop equations. 
After some straightforward algebra, we find 



i^>(t,x;i',x') = (Ob|^(i,x)V;(t',x')|Ob) = 



(2^) 
(jl+m) cos (3p 



iE„t' 



smfip e'^^j5lo Spe 



Upon reordering of the terms we find the Green's function quoted in ( IV.3| ). 

A similar calculation leads to the transformed Green function S^{t, x; t' , x') which is defined as 

-z5<(t,x;t',x') = (Ob|V^(t',x')^(t,x)|Ob) , 



Following the same steps leading to ( [A.3| ) we find 

-z5<(t,x;i',x') = (Ob|^(i',x')^(t,x)|00 = 



sin/3p e 



-i5rj ^iEjjt 



Sp + COS f3p 7570 e 



-iE„t' 



which upon reordering of terms gives the form quoted in expression ([V.4) 



(A.3) 



APPENDIX B: BOGOLIUBOV TRANSFORMATION AND TADPOLE DIAGRAM FOR 

INHOMOGENOUS SYSTEMS 



In this Appendix we generalize the Bogoliubov transformations described in the homogeneous case to the case of 
inhomogeneous condensates. Unlike the homogeneous case in which the generator of the Bogoliubov transformation 
creates particle-antiparticle pairs of zero total momentum, in the inhomogenous case the total momentum of the pair 
is non-zero. 

Consistent with perturbation theory we now find the corresponding Bogoliubov transformation to lowest order in 
the Yukawa couphng, thus cos Pps « 1 , sin/3ps « (3ps — 0{g). Since the self-energy is already of 0{g'^), to lowest 
order we only need to focus on the tadpole term J[,(x, t). 

The Bogoliubov transformation in lowest order reads 



/cPp' 

dH-p.s) = - j ^_||_p:,^(p',p)6,(p',.') + 4(-P-^) 
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with pss' ~ 0{g). With this choice the transformation leaves the canonical anticommutation relations unchanged up 
to terms of order p^^, . In order to compute the transformed Green functions we need the expectation values of bilinear 
combinations of creation and annihilation operator. We find the following expectation values that are necessary to 
compute the Green's functions , 



(0b|6(p,,sK-p',s')|0b) 



Pss' (P, P') • 



Ps"s'{p",p')] {0^Mp,s)bl{p",s"M)d,s'{27rf2EpS^{p-p") 



(0^|6(p,s)6t(p',s')|0^) = {27Tf2EpSs,'S^ip-p') , 
(Ob|dt(-p,s)d(-p',s')|0,) = 0, 
(0,Mt(-p,s)6t(p',s')|0,) = -p:,,(p',p) . 

(Ot|d(p',s')dt(p,s)|Ot) = {2TTf2Eo6ss'S''{p-p') , 
(Ob|6t(-p',s')6(-p,s)|Ob)=0, 

W{-p',s')b{p,s)\0,) = pss'{p,p') , 

{0,\b\p',s')d\-p,s)\0,) = p:,,(p',p) . 



This yields the transformed Green's functions, 

i5>(t,x;i',x') = (Ob|V(i,x)Vi(t',x')|Ob) = 

d^p 



+P..'{P,P'} f (P,») V{-p',,') e'f^^.'+i^.''')} . 



and 



-i5<(i,x;i',x') = (Ob|Vi(t',x')V(i,x)|Ob) 



/ / {p:.(p'.p) n-P,«) ^(P',^0 e-(-^*-^'*') 

+ Pss' (P, p') U{p, s) V{-p', s') e^(^.*+^.'*') } . 
to first order in p.,,,/ (p, p'). 

Using these results we now can evaluate the tadpole graph in the inhomogenous condensate. That is, the expectation 
value of < ■^V > which plays the role of an external current in the equation of motion. Inserting the above explicit 
expressions we find 

Jt(t,x) =*5tr5>(t,x;<,x) =5;^ y tr [C/(p, s)[7(p, s)] 

_ V / / _fp_ i{p-p')^ 

j^J {27^f2Ep J {2T,f2Ep, 

{p:,,(p',p) tr [V{-p,s) C7(p',.')] e^(^''+^^')* 
Pss'M) tr [U{p,s) V{-p',s')\ e-(^''+^=P')t} . 
The traces over the spinors yield. 
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tr V{-p, s)U{p\ s') = U{p', s')V{-p, s) 



(B.l) 



^ {Ep ~ m){Ep, +m) + ^{Ep + m){Ep, - m) 



tr C/(p, s)F(-p', s') = Vi-p', ,s')C/(p, s) 



^J {Ep - m){Ep, +m) + ^{Ep + m){Ep, - m) 



We see that the relevant part of pss'(p,p') contributing to eq.(BJ) is odd in s and diagonal in ss'. We therefore 
choose, 



Pss'{p,p') = 7(P,P') 



Sss' 



y^{Ep-m){Ep, +m) + ^{Ep + m){Ep, - m) 



Then, 



{2tt)^2Ep 



Jb(t,x) = Amg 



{2nf2Ep 



-i{Ep+E. 



The first term is again space and time independent and is absorbed into a shift of the condensate. The second term 
displays space and time dependence in a factored form. It will be used to compensate for the initial time singularities 
of the self-energy. 

APPENDIX C: ANALYSIS OF THE SELF ENERGY KERNEL 

In this Appendix we provide the details for the various contributions to the self-energy. The integrals that enter in 
the expression for the self energy kernel can be related to the following one defined in dimensional regularization 



1 



E++E_ 



(2^)3- 2E+E^ {E++E-f^ql 
d^-'p 1 



i2nr 



{p — q/2)^ — +10 {p + q/2)^ — -f io 



167r2 



da In 



m 



ni^ + a{l — a) (q^ — q^) 



Where we have introduced the shifted momenta p± = p ± q/2 and energies E± = y pj_ + m?. is defined as 

Le = 7 + In 5- • 

We now consider the various integrals defined in section VIII. In doing so we will shift the integration variable p so 
that p p+ = p -I- q/2 and p' = p — q ^ P- = P — q/2. Then the numerator arising from the Dirac trace takes the 
form 



E+E_ + p+p_ -w? ^ Ue+ + E.f -2[m^ + \- 



We then have 



Si(q^)- -85^/1 (q2) 



with 
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d? E+E_ + p+p_ - 
(27r)3-^ iE+E^ {E+ + E^) 



1 



1 



8 J (27r)3-«^ 

The first integral, including the prefactor, is equal to 

(f-^p 1 



4 



1 



(27r)3-^ 2£;+i;_(jE;+ + £;_) 



r(ie + l) 



(27r)3~^ AE 327r2 
The second integral is the basic integral I{q^,c^) at = 0. Altogether we obtain 

Ei(q') = Ei(0) + q's;(0) + AEi(q2) 

with 



Si(0) 



3rn 



2^2 



47r2 



S;(0) = -<5Z=|-^L, 



AI]i(q^)=(m^ + :V)|^ 



da In 



vn? + a(l — a) q^ 



Here we have introduced the renormalization constants corresponding to a renormalization aX, — Q. 
For S3(q^) we have 



with 



/3(q') = 



d3- 



S3(q') = 8 5'/3(q') 

1 (i;++i;_)V2-2(m2 + q2/4) 



(27r)3-M£;+£;_ {E++E.Y 
This integral can be related to the integral /((?o: Q^) ^'^d its derivative w.r.t. q^, at go = 0. We find 

S3(q') = -<5Z + AE3(q2) 



where (5Z has been defined in eq.C.l). The finite part is 



^27r2 



da In 



1 + a(l — a) — 2 



da- 



rn^ 
a(\ — a) 



m2 + a{l — a)q2 



(C.l) 



From the way in which we have introduced 5Z in Si and S3 it is apparent that the covariant counterterms 5Z[(j)+c^ 
in the equation of motion will absorb these divergences. 

We finally consider the Laplace transform of the subtracted self energy kernel introduced in eq.( VII.l ) 

1 (£;+ + £;_)2/2-2(m2 + q2/4) 



E^ {E+ + E^f[{E+ + E^Y + s'^] 



Comparing with the standard integral /(qq, q^) we see that besides the continuation to the Euclidean region, q^ —s^ 
we have additional denominators. These can be obtained via subtraction. Wc have 



d^p 



(27r)3 2E+E^ {E+ + E^)[{E+ + S_)2 + s^] 



= --^[/(-.^q^)-/(0,q^)] 



1 1 
^167r2 



da In 



a{l — a)q-^ 



m2 + q;(1 — Q;)(q2 + s^) 
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We proceed analogously for the second integral and finally obtain 

„2 (-1 



27r2 



da 







1 



a{l — a) 
4 J m? + a{l — a)q2 

The a integrations can be performed analytically. 



In 



rr? + ol(\ — a)c^ 
m? + ol(\ — a)(q2 + s^) 



(C.2) 



APPENDIX D: DETAILS OF THE ANALYTIC SOLUTION 

We have obtained in section the solution of the equation of motion and its solution via Laplace transform. We 
consider at first the unrenormalized equation. The solution reads 



V'q(s) 



S ,^q(0) + 0q(O) 



so that 



2-Ki 



-ioo+c 



s2 + Af2+q2+a(s2,q2) ' 

S0q(O) + (/)q(0) 



— lOO + C 



s2 + AP + q2 + CT(s2,q2) 



As usual [ [I5yi6| we shift the contour to the left so that finally it includes the cuts, and eventually poles, on the 
imaginary s axis and a circle at \s\ oo around the left half, which does not contribute for positive t as the 
exponential exp st tends to zero there. In doing so we make use of the causality condition that there are no zeros in 
the left half of the complex s plane, as required by causality. Along the cut s^t s = iuj with l^m^ + q2 < < oo we 
define the real and imaginary parts of the kernel a by the convention 

a((iw±e)2,q2) ^a,^(-c^^q2)±ia/(-w2^q2) . 

As a only depends on this also fixes the relative signs of the imaginary parts of a on the lower cut for which 
— oo < w < —2^J'm? + q2. We then obtain 



'-^'^ - 2^. 



i duj e"^' disc ■ 



1 

'27ri 



disc 



-^2 + Af2 + p2 + ct(-w2 ± ie, q2 

-za;</.q(O) + 0q(O) 
-cj2 + M2 + p2 + a(-w2 =F le, q2) 



with LOc = 2 •\/q2 + w? for the two fermion cut. The spectral density is obtained from the discontinuity across the cut 

1 



S{ijj, q) = i disc 



-cj2 + M2 + p2 + a(-cj2 + ie, q2 



2aj(-c.2,q2) 



M2 



+ o-fi(-ti;2 + ie, q2)] + a'j[-uo'^ + ie,c^) 



(D.l) 



In the case in which the scalar particle is unstable, i.e. M > 2m there is a resonance above the fermion-antifermion 
threshold and no support for the spectral density below threshold However in the case M < 2m the scalar 

is stable and cannot decay, now the spectral density has support above and below threshold. Below threshold the 
spectral density is a delta function at the position of the renormalized pole, to include the stable pole below the 
two particle threshold in the description we now define ujc = 0"^ to distinguish that the origin is excluded from the 
integration region. The pole in the stable case is obtained from the identity 



nSi-LO^ + + p2 + dRi-LO^ + ze, q^)) 



so that 
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duj cos ujt 



0+ 



[-uo'^ + A-P + p2 + ani-uj'^ + ie, <f)] + ct|(-cj2 + ie, q2) 

0q(O) a/(-Cc>^q^) 

Af2 + p2 + CT/?(-u;2 + ie, + aj{~uj'^ + ie, q2) ' 



(D.2) 



In order that this equation and its time derivative be consistent at i = we have to require the sum rule 



duj 



uj a/(-w2,q2) 



Jo+ 



[-cj2 + M2 + p2 + afl(-a;2 + ie, q^)]' + dji^uj^ + ie, q 



In order to derive this sume rule we require, as already mentioned above, that the denominator s'^ + q2 + AP + a{s'^,q^) 
has no zeros in the left half of the complex plane. We have to assume furthermore that (t(s2, q2) increases less strongly 



as s as \s\ — > oo in the left half of the complex plane. Under these assumptions we have the identity 

J_ /■ . s 

in 



q2 + M2 + CT(s2^q2; 



= 



if the integral is carried out along the contour enclosing the left half of the complex plane. The contour consists of 
an integral along the left of the imaginary s axis and a semicircle at \s\ = oo. The latter one contibutes 



1 

in 



ds 



= -1 



The integral along the imaginary axis is given by 



duj 



duj 



- q2 + A/2 + 5-(-tj2 _ jo) 

;2 + 1/2 + p2 + ct(-w2 _|_ q2 



s2 + q2 + Af2 + CT(-a;2 + jo) 



(D.3) 



The two parts of the contour integral have to add up to zero, which yields the sum rule. 
For the renormalized equation of motion we rewrite the result obtained in section VII as 



ioo+c 



ds 



ioQ-\-c 



s0q(O) + 0q(O)] Fi(s2,q2)+0q(O)iF2(s2,q2) 



with 



^q(^) = 7^ 

2m 



The function Fi{s'^, q2) has analyticity properties analogous to the fraction l/(s^ + q^ + A/2 + a) considered above. 
In particular the discontinuities along the positive and negative imaginary axis have the same relative signs, it has 
no singularities in the left half of the complex plane, and the limiting behavior as |s| oo is For the first 

property it is essential to note that a-g only depends on the square of the variable s, see (C.2). For the last property is 
sufficient to notice that as{s'^,q^) behaves as lns2 as \s\ —> oo, so the terms proportional to as dominate in numerator 
and denominator. The function F2(s2,q2) has analyticity properties analogous to those of Fi{s'^,q^), and decreases 
asymptotically as l/s2. The relative signs of the imaginary parts along the cuts arc the same as for Fi(s2^q2) 
Collecting all terms wc find 





H 


- AS3(q') H 


h s 


a(,s2,q2) 




S2 


"l + AS3(q2)- 


hCT,(s2,q2) 


+ 


q2 + A/2 - 


HASi(s2,q2) 






CTs(s^q 








S2 


"l + AS3(q2)- 


h CTs(s2,q2) 


+ 


q2 + A/2 - 


hAEi(q2) ■ 



2 f°° 

f, it) — — j dcj Tcosuji 0q(O) w Im/^i(— cj^ + io, q^) 



(D.4) 



-sinujt 0q(O) ImFi(- 



io,q^) 



,(0) 



cosLut ImF2(— ti)^ + io, q^) 
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For the consistency of the left and right hand sides and their first derivatives w.r.t t the sum rule 

2 f°° 

1^— duj uj linFi{-uj'^ + io,q^) 
has to be satisfied, ft follows again by considering the integral 

— idss Fi{s^,c^) 

ITT J 

along a closed contour around the left complex half plane. One needs, furthermore, 



(D.5) 



= -/ dcj - ImFsf-cj^ + io,q^) 

which follows from analagous considerations, using that in this case the infinite semicircle does not contribute as the 
integrand behaves as there. 

We next have to consider the term proportional to <?!'q(0) and the second time derivative of eq.( D.4 ). From the 
renormalized equation of motion (VI. 2) at t = we derive immediately 



?iq(0) = -0q(O) 



q2 + Af2 + AEi(q2) 
l + AE3(q)2 



The second derivative of (D.4) at t = reads 

2 



i(0) = - 



duj 



0q(O) uj^ liTvFii-ijj^ + io,c^) + u? 0q(O) Im F2 (-tj^ + io, q^) 



We express, on the right hand side, ^q(O) by ^q(O), using ( |D.6D . Then we can write the sum rule as 



1 = - 

TT 



duj LO Im 



l + AE3(q2) 



92 + M2 + AEi(q2) 
The expression in the square brackets can be written explicitly as 

(1 + AS3 + CT.) 1 + AEi 



+ i^2(-w^ + jo,q^) 



-uj^{l + AS3 + + q2 + Af2 + AEi + M"^ + ASi 



(D.6) 



(D.7) 



^cj2 \l + AS3 + CT, ) + q2 + M2 + ASi 
1 + ASi 



q2 + M2 + AEi 



1 + AS3 + a. 



-tj2 1 + AS3 + a J + q2 + AP + ASi 



Only the imaginary part if this expression occurs in the in tegr and. So the first term on the right han d side does not 
contribute, and the second term is just Fi. The sum rule (D.7) reduces therefore to the first one(|D.5|). 
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FIG. 2. Stable case: M = 1 ; m = 
and normalized so that J d^x4>{0, x) 



1 ; g 
= 1. 



1 with a gaussian profile for the condensate at t = 0, given by ( VII.2| ) with Ro — 1 
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